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Extended phase space (EPS) formulation of quantum statistical mechanics treats the 
ordinary phase space coordinates on the same footing and thereby permits the defi- 
nite the canonical momenta conjugate to these coordinates . The extended lagrangian 
and extended hamiltonian are defined in EPS by the same procedure as one does for 
ordinary lagrangian and hamiltonian. The combination of ordinary phase space and 
their conjugate momenta exhibits the evolution of particles and their mirror images 
together. The resultant evolution equation in EPS for a damped harmonic oscillator, 
is such that the energy dissipated by the actual oscillator is absorbed in the same rate 
by the image oscillator leaving the whole system as a conservative system. 
We use the EPS formalism to obtain the dual hamiltonian of a damped harmonic os- 
cillator, first proposed by Batemann, by a simple extended canonical transformations 
in the extended phase space. The extended canonical transformations are capable of 
converting the damped system of actual and image oscillators to an undamped one, 
and transform the evolution equation into a simple form. The resultant equation is 
solved and the eigenvalues and eigenfunctions for damped oscillator and its mirror 
image are obtained. The results are in agreement with those obtained by Bateman. 
At last, the uncertainty relation are examined for above system. 



1 Introduction 

Although, the formulation of dissipative systems from the first principles are cumbersome and 
little transparent, however, it is not so difficult to account for dissipative forces in classical 
mechanics in a phenomenological manner. Stokes' linear frictional force proportional to the 
velocity v, coulomb's friction ~ ^, Dirac's radiation damping ~ v and the viscous force ~ V^v 
are noteworthy examples in this respect. Unfortunately, the situation in much more compli- 
cated in quantum level (see Dekker,1981, [1], and the references there in ). In this review article 
on classical and quantum mechanics of the damped harmonic oscillator, Dekker outlines that: 
"Although completeness is certainly not claimed, it is felt that the present text covers a sub- 
stantial portion of the relevant work done during the last half century. All models agree on the 
classical dynamics ... howover, the actual quantum mechanics of the various models reveals a 
considerable variety in fluctuation behavior. ... close inspection the further shows that none of 
them ... are completely satisfactory in all respects. "As an example of the dissipative systems, 
the damped harmonic oscillators is investigated through different approaches by different peo- 
ple. Caldirola [2] and Kanai[3] using the familiar canonical quantization procedure, obtained 
the Shrodinger equation which gives the eigenvalue and eignfunctions for damped oscillator. 
However the difficulty with this approach is that, it violates the Hisenberg uncertainty relation 
in the long time limit. Another approach is the Shrodinger- Langvin method which introduces 
a nonlinear wave equation for the evolution of the damped oscillator [4]. In this method the 
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superposition principle is obviously violated. Using the Winger equation, Dodonov and Manko 
introduced the loss energy [5]. As consequence of the dissipative Bateman, by introducing a dual 
hamiltonian considered the evolution of the damped oscillator in parallel with it mirror image 
[6] . The energy dissipated by the actual oscillator of interest is absorbed at the same rate by the 
image oscillator. The image oscillator, in fact, plays the role of the physical reservoir. Therefore, 
the energy of the total system, as a closed one, is a constant of motion. 

Here we use the EPS method [7] to investigate the evolution a damped harmonic oscillator. 
The method looks like the Bateman approach, however, the uncertainty principle, when looked 
upon from a different point of view, is not violated. That is, the extended uncertainty relation is 
satisfied for combination of actual and image oscillators, while reducing into ordinary uncertainty 
relations for actual and image oscillators in zero dissipation constant limit. 

This paper organized as follows. In section 5, a review of the EPS formulation is given. 
In section 7, we investigate the quantization procedure for the damped harmonic oscillator. In 
section 4, we use the path integral technique directly to calculate the exact propagators, and then 
the uncertainties of position and canonical momenta for the actual and mirror image oscillator 
system. Section 5 is devoted ko concluding remarks. 



2 A review of the EPS formulation 

A direct approach to quantum statistical mechanics is proposrd by Sobouti and Nasiri[7], by 
extending the conveltional phase space and applying the canonical quantization procedure to 

extended quantities in this space. Assuming the phase space coordinates q and p to be indepen- 
dent variables on the virtual trajectories, allows one to defined momenta vr^ and tt^, conjugate 
to q and p, respectively. This is done by introducing the extended lagrangian 

^<i,p,q,'p) = -qp-pq + ^''{q,q) + ^'''{p,p) (i) 

where and are the q and p space lagrangians of the given system. Using Eq. (1) one may 
define the momenta, conjugate to q and p, respectively, as follow 

P, (2) 

- q- (3) 



dC 






dq 


dC 






dp 



In the EPS defined by the set of variables {q,p, tt^, tt^} , one may define the extended hamiltonian 

Ti-{q,P,'Kq,'ITp) = qTTq+pTTp- C = H{p + TTq,q) - H{p,q + Trp) 

E}_c&^^ t _ 
n^Opn^i dq^ ^' 

where H(q,p) is the hamiltonian of the system. Using the canonical quantization rule, the 
following postulates are outlined: 

a) Let q, p, TTg and tt^ be operators in Hilbert space , X, of all square integrable complex 
functions, satisfying the following commutation relations 

d 

[iTp,p]=-ih, iTp = -ih-^, (6) 

= [7rg,7rp] = 0. (7) 

By virtue of Eq. (5), Eq. (6) and Eq. (7), thhe extended hamiltonian, H, will be an operator 
in X. 



S. Nasiri and H. Safari 



3 



b) A state function xiQ^P^t) G X is assumed to satisfy the following dynamical equation 

dx d d 

^^~dt ^ '^X=[H{p-ih—,q)-H{p,q-ih—)]x 

The general solution for this equation is 

x{q,P,t)=m^*ip)e-^''', (9) 

where ipiq) and (f)(j)) are the solutions of the Schrodinger equation in q and p space, respectively, 
c) the averaging rule for an observable 0{q,p), a c-number operator in thij formalism, is given 
as 

<0{q,p)>= J 0{q, p)x* {q, P, t)dpdq. (10) 
for details of selection procedure of the admissible state functions, see Sobouti and Nasiri[7]. 

3 Damped harmonic oscillator in EPS 

Extended hamiltonian of Eq. (4) for undamped harmonic oscillator is given by 

^ = ^^9^ + P'^q - \'^P^ - IT^P- (11) 

By a canonical transformation of the form 

Q = Ql,'^q = -TTgi -Pi, 
p = pi,Trp = -TTpg - qs, 

Equation (10) yields 

n = lnl+qf-lnl-lpl (12) 

This extended hamiltonian evidently represents the subtraction of hamiltonians of two indepen- 
dent ideotical oscillators, which is called actual and image oscillators [5]. The position q and 
momentum Tig denote the actval oscillator while, p and TVp denote the image oscillator. The mi- 
nus sign has its origin from Eq. (4) and has an important role in this theory [7]. The following 
canonical transformation 

q2=jl, ■n-qi = TTg^ + Ml, 
P2 =Pl, TTpi = TTpj - Xpi. 

(13) 

changes the extended hamiltonian of an undamped harmonic oscillator into that of the damped 
one, i. e. 

^2 = ^{tt^ + 2Ag27r<,2 + oj^qh - \{t^I^ - ^XpoTTp, + co^pI}, (14) 
where co = 4 + jX. One further transformation generated by 

^2(92,P2,7rg3,7rp3) = g27rg2e~'^'* +P27rp3e'^* (15) 
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finally leads to 

^3 = 21^53^ +^Q e f-^i^PB^^ +UJ p3e |. (16) 

The first part of the extended hamiltonian Eq. (11) is Caldirola- Kanai hamiltonian, which is 
widely used to study the dissipation in quantum mechanics [3]. Using Eq. (15), the extended 
Hamilton equations [7] gives the following evolution equations for actual and image oscillators, 
respectively 

q's + OAga + cj^Qs = 0, (17) 

and 

PS - 2Xp3 + J^ps = 0. (18) 

Almost trivially, the energy dissipated by actual osciluator, with phase space coordinates (^4, -Kq^) 
is completely absorbed at the same pace by the image oscillator with phase space coordinates 

(P0,vrp3). 

As usual, the dynamical variables (^3, tt^^) and (ps, VTpg) are considered as operator in a linear 
space. They obey the commutation relations Eqs. (5) - (7). The dynamical equation, Eq. (8), 
now becomes 



S 2^-2Xt I , ,0 2^2At\ ^ r„ 2 2\t , , ,2 2^-2Atx\,, n Q^ 



2-2Xt , , ,0 2^2At\ 1 

By an infinitesimal canonical transformation which in quantum level corresponds to the following 
unitary transformation 

U = exp(^{e'^*92' + e-^^'p2} + '-^{qsTTg^ - PsTTp^}), (20) 

equation (10) may be written as 

dx 

ih^ = nx 

= {ns = ^{7r,3^ + u;'\s'} - l{7r,J + u'W}.}x. (21) 
where w' = w + i\. Eigenvalues of Eq. (22) are [7] 
f — E — E 

= (n - m)huj'. (22) 
Which are in agreement whit those obtained by Bateman. The eigenfunctions of are 

X{q?„p?„t) = ^{q^W{pz)e-^ni-^P\ (23) 

where ijj{q) and 0(p) eigekfunctions of harmonic oscillators in configuration and mymentum 
space (Hermit functions). Then the eigenfunction for Eq. (12) reads 

x' = Ux(.q3,P3,t) 

= expi^ie^^^qs'^ + e-^^W})expC-^{q3Trg^ - P3Trp^})x 

= expC^{e^^%^ + e-^^W}me^'q3)4>*{e-^%)e'iP''iK (24) 

where , are unitary transformation of two consecztive canonical transformation. Finally, eigen- 
function of Eq. (15) reads The above eigenfunction are completed to investigate the uncertainties 
relations for the combined systems actual and images oscillators system. 
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4 Uncertainty relations for actual and image oscillators 

In this section we calculate the uncertainties in position and momentum for the actual and the 
image oscillators. We calculate the extended propagator [9] for the combined actual and the 
image oscillators as follows 

K{q,p,t,qi,pi,ti) = 



'dmh sinvo'it — ti) 

2 smLO [t — ti) ijj' 

+ e~^^^'^'^qi{cosuj'{t - U) + ^sin{io'{t - U)) - 2qqi}] 

A, u;'e~'^(*+*') ^ r \(t t\ I, s ^ , I, 

X exp - - — X {e"^(*"*»V cosw t - U) + —sin{J(t - U)) 

2 smLo'[t — ti) to' 

+ e^(*-*')pf (cosw'(t - t^) - —sin{J{t - ti)) - 2ppi}]. (25) 

When A — > then Eq. (|25jl reduces to the familiar form of the undamped extended harmonic 
oscillator propagtor [8] . We assume that the initial state function for combine system in ground 
state is Xooil^P^ 0) = (7rJ^)~2exp(— '^^"^f ) where delta is the width of the extended wave packet. 
Then one gets using Eq. () 

vr 1 ioj' cosuj't ^ 

7T)^e.p[-^{-,e2^*(l + [-(-)2 + 2(-)^-l] 



UJ UJ 

sm UJ H sin2uj t) — ^i^: x l\cosuj t smuj t 

uo' h sinuj't uj' 



I A ,,,, l,?i,n ,A,9 , , A i Lo' cosuj't — 

( cosJt + -sinJt)][l + - - 2 + 2-2 _ l]sin2Jt + - sin2J t]]-^][-- + --{—— + - 
uj' 0* uj' ijj ijj 2d^ 2 h smuj't u; 



9 / — ^ ■ rt I \— 1 -r^ ^ /r / ^ • / 

s^n w H -sin2u!t) +«|^: — x ([costj t H -smujt 

A . /xTr- 1/^n9 ^/A,9 ^t.^/ a 



( cosw't -sinuj't)][l + ^{^Y + 2(— )2 - l]sin2Jt -sin2uj't}}-^\e-— . 

uj' 0* Uj' u>' u> 

Using Eq. Q the uncertainties of positions and momenta can be computed for the actual and 
the image oscillators as follows 



< Aq> = V< g2 > _ < q >s 



^e-''{l + [{^)\^f + (4)^ - l]«mVi + ^s^n2u;'tf^. (27) 



and 



< AvTg > = J< 7r2 > - < 7r„ >2 



Ae-A*{i + [(^)4(1)2 + (A)2 _ _ ^s^n2u:'t}h. (28) 
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and 

< A^? > = \/< ■p^> - <P>'^ 

= ^e^Hl + [C^)\-f + {-,? - l]sM - ^sin2u't}h. (29) 

and 

= 4=e'Hl + + (A)' - Mn^^'t + A,i„2u;'t}i, (30) 

The above results for actual and image oscillator, in separate form, is in agreement with those 
obtained by Bateman. However, when A 7^ 0, it is not possible to separate the oscillators 
as Bateman docs and the Heisenberg uncertainty relations would not hold for each oscillator, 
separately, (see figs. 1 and 2). In fact the uncertaintyrelations would also looked upon for 
combined system n extended phase space as shown in fig. 3. 
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